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We investigate the effects of including the two-loop QCD correction to the static short distance 
1/r potential on the spectroscopy and decays in the charmonium and upsilon systems. We use a vari- 
ational technique with the full two-loop corrected potential to determine a set of unperturbed trial 
wave functions and treat the relativistic and one-loop corrections as perturbations. This approach 
results in more accurate descriptions of the hyperfine splittings in both the 66 and cc systems. 
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In an earlier publication [lj, two of us investigated the spectrum, leptonic widths, and electromagnetic decays 
of the charmonium and upsilon systems. We adopted a variational approach based on determining a set of trial 
wave functions derived from an unperturbed Hamiltonian consisting of a relativistic kinetic energy term, a long 
43 " range confining potential and a short distance 1/r potential that included the one- loop QCD correction. These 
wave functions were then used to compute the contributions of the remaining relativistic and one-loop corrections 
^ • to the short distance potential using perturbation theory. The resulting masses of the various states were then 
QJ , compared with a subset of the experimental values using a x 2 test and the procedure was repeated until a % 2 
43 minimum was obtained. 

The model in Ref. [IJ yielded quite reasonable results, but there were some areas of concern. These had to do 
with the size of the hyperfine splittings and the values of some leptonic widths. Since both of these observables 
are sensitive to the values of the radial wave functions at the origin, we suspected that the trial functions, even 
with the inclusion of the one-loop QCD correction, were not sufficiently well determined at short distances. 
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CO 1 To address this problem, in this paper we include both the one-loop and two-loop corrections to the short 

^sO \ distance 1/r potential in the unperturbed Hamiltonian. For the full static QCD potential to two-loops, we use 
the results of Ref ■ The two- loop derivation was performed using the MS renormalization scheme. As discussed 
below, the use of the MS scheme requires a slight modification to our one-loop perturbative potential. With 
CNI ■ this modification, we can then repeat the variational calculation previously mentioned. 

In Sec.[H we discuss the modifications of the perturbative potential associated with changing it to the MS 
renormalization scheme and, in Sec. [31 we present the results for the charmonium and upsilon systems. Our 
, conclusions are contained in Sec.|4j Several calculational details are presented in the Appendices. 

2. MODIFIED SEMIRELATIVISTIC MODEL 

We use a semi-relativistic unperturbed Hamiltonian of the form 



H = 2m v / p 2 +m 2 + Ar + V(r) , (1) 
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where V(r) is the two- loop QCD short-range static potential Q. It is calculated in full detail in Appendix A 
and is 
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where \i is the renormalization scale and rif denotes the number of light quarks. The perturbative potential has 
the form 



H' = V L + V S . 

Here, Vl is the order v 2 /c 2 correction to the confining potential, 
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The short distance potential Vs 



V s = V HF + V LS + V t + Vsi 



(3) 
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must be modified in order to be consistent with the two-loop QCD static potential V(r). This amounts to 
changing the renormalization scheme from the Gupta-Radford version used in [l| to the MS version used in 
obtaining V(r). The transformation between the as used in [l[ and as is given by 



_ a s 49 10 
a s = a s [l + —(— ~ —n f )} . 

From this relationship we now can write down the modified short distance potential. 
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3. RESULTS 



Although we are treating Vl and Vs perturbatively, we will retain the (5-functions in Vrf and Vsi rather 
than taking the usual approach of 'softening' them by using a sharply peaked, differcntiablc function. This is a 
departure from [l[. However, with the modified Hq, we found that the softened 5-function offered no substantial 
improvement to the overall cc and bb spectra, and had the disadvantage that it predicts a considerably smaller 
T(1S) — 77b (IS) hyperfine splitting than the one observed 0-Q along with a less successful description of the 
leptonic widths. The final parameters that resulted from the variational calculation are shown in Table [IJ 
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cc 


bb 


A (GeV 2 ) 


0.149 


0.168 


as 


0.304 


0.213 


m q (GeV) 


1.69 


5.03 


At (GeV) 


1.95 


8.36 



TABLE I: Fitted Parameters for the cc and bb systems 

3.1. Charmonium System 

The results for the spectrum of the charmonium system are shown in Table [II] We set fij to 4 for the 
calculations. The description of the cc spectrum with the modified non-perturbative wave functions and the 



m c e (MeV) 


Theory 


Expt 




2978.7 


2980.4 ± 1.2 


j/4>(isy 


3095.5 


3096.916 ±0.011 


xMipy 


3419.4 


3414.74 ± 0.35 


Xel(lP)* 


3505.9 


3510.77 ±0.07 


Xe 2 (LP)* 


3556.2 


3556.20 ±0.09 


h c (lP) 


3525.9 


3525.93 ±0.27 


Vc(2S) 


3603.2 


3638.0 ±4.0 


V(2S)* 


3690.9 


3686.093 ± 0.034 


1>(1D) 


3802.2 


3771.1 ±2.4 


1 3 D 2 


3821.4 




1 3 D 3 


3831.0 




1 X D 2 


3822.5 




Xco(2P) 


3843.6 




Xci(2P) 


3928.1 




Xc2(2P) 


3978.9 


3929.0 ± 5.4 


h c (2P) 


3948.6 




r? c (3S) 


4006.9 




^(35) 


4087.8 


4039.0 ± 1.0 


^{2DY 


4148.0 


4153.0 ±3.0 


2 A D 2 


4170.5 




2 3 D 3 


4184.0 




2 1 Z? 2 


4172.8 





TABLE II: Results for the cc spectrum are shown. The states denoted by a * are used in the fitting pro cedure. In order 
to account higher order corrections, we used a 2.5 MeV error offset in the fit. All data is taken from [lGl ]. 

retention of the 5-functions in the perturbation is consistent with our earlier calculations [l|. The result of 116.7 
MeV for J /ip{\S) — r/c({S) splitting compares favorably with the 116.6 ± 1.2 MeV experimental value. The 
corresponding iJj(2S) — rjc(2S) result of 87.6 MeV remains significantly larger than the experimental value of 
48.1 ± 4.0 MeV. On the other hand, the two- loop unperturbed wave functions provide a much better description 
of the leptonic widths, shown in Tablc lTOl The dipole decay rates, which depend on the matrix elements of r, 
are not particularly affected by the changes in the short distance behavior associated with the inclusion of the 
two-loop corrections. 
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i e e ^Kev ) ineory 


ILXpt 




(J.UU ZC U. ±"4: 


^(25) 2.80 


2.48 ± 0.06 


i/^S 1 ) 2.00 


0.86 ±0.07 


V»(l-D) 0.14 


0.242 ± 0.030 



TABLE III: Leptonic widths of the ip(nS) and the ip(lD) states are shown. The ip(nS) widths include the QCD correction 
factor (1 — 167tqs/3). 



3.2. Upsilon System 



The results for the 66 spectrum, computed using rif = 5, arc shown in Table fVl For the levels that are not 
particularly sensitive to the short distance behavior wave functions, the results are compatible with Ref. [l[. 
The main improvement is seen in the T(IS') — r]b(lS) hyperfme interval fllT - fl3T ]. which has increased from 47.0 
McV to 67.4 MeV, much closer to experimental value of 69.5 ± 3.2 MeV. The 2S and 3S hyperfme intervals are 
also larger than those in [l[. Our calculation is also consistent with the recent discovery of evidence for a 3 3 Pj 
multiple! from the ATLAS detector (IHl^. The leptonic widths are given in Table [TV] 



V ee (keV) 


Theory 


Expt 


T(1S) 


1.340 


1.340 ±0.018 


T(2S) 


0.644 


0.612 ±0.011 


T(3S) 


0.486 


0.443 ± 0.008 


T(4S) 


0.444 


0.272 ± 0.029 



TABLE IV: Leptonic widths of the T(nS) states are shown. The tjj(nS) widths include the QCD correction factor 
(1 - 167ra s /3) 



4. CONCLUSIONS 



The inclusion of the two-loop corrections to the short distance static potential used in the determination of the 
unperturbed trial wave functions improves description of cc and 66 observables that are sensitive to the behavior 
of the wave function near r — > 0. This enabled us to retain the (5-function terms in the perturbativc potential 
and account for the larger ground state hyperfme splittings in both the charm and upsilon systems. These 
improvements can be obtained without a dramatic change in the descriptions of the overall cc and 66 spectra. 



Appendix A: Two- loop short distance potential calculational details 



The momentum space version of the two-loop short distance potential, V{k 2 ), in terms of as is given in the 
paper by Y. Schroder 0] using the MS renormalization scheme. In order to calculate the coordinate space 
potential, V(r), it is necessary to evaluate integrals of the form 
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m bb (MeVJ 


Theory 


r^xpt 


VbO-S) 


9391.5 


9390.8 ± 3.2 


T(15)* 


9459.0 


9460.30 ± 0.26 


X60(lf)* 


9865.7 


9859.44 ± 0.52 


X6l(lP)* 


9892.7 


9892.78 ± 0.40 


Xb2(lP)* 


9910.9 


9912.21 ± 0.40 


h b (lP) 


9900.2 






9984.6 




T(2S)* 


10021.3 


10023.26 ± 0.31 


T(1D) 


10151.1 




1 3 D2 


10157.2 


10161.1 ± 1.7 


1 3 L>3 


10161.5 




1 D 2 


10158.1 






10231.3 


10232.5 ± 0.6 




10253.8 


10255.46 ± 0.55 


xt2{2py 


10269.1 


10268.65 ± 0.55 


h b (2P) 


10260.2 




1 3 F 2 


10352.9 




1 3 -F3 


10354.7 




1 3 F 4 


10355.2 




1 F 3 


10354.5 




t? 6 (3S) 


10331.8 




T(3S)* 


10361.2 


10355.2 ± 0.5 


Y(2D) 


10442.2 






1 (Til 3 7 




X6i(3P) 


10534.4 




Xf> 2 (3P) 


10548.7 




(3 3 Pj) 


10540.0 


10530 ± 10.3 



TABLE V: Results for the bb spectrum are shown. The states denoted by a * are used in the fitting procedure. In order 
to account higher order corrections, we used a 2.5 MeV error offset in the fit. All data taken from [Icfl . The notation 
(3 3 Pj) denotes the spin average of the 3 3 Pj levels. 



The latter integral is readily evaluated using Mathematica for the values n = 0, 1, 2 that are needed to obtain 
V(r). Explicitly, the expression for V(k 2 ) is Q 
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where c\(x) and 02(2:) are given by [l9j . 

ci(x) = cii - p \n(x) (A4) 
c 2 (x) = a 2 +/3 2 ln 2 (x)-(^i+2/3oai)ln(.T). (A5) 

The parameters a±, a 2 , /?o an d Pi arc given in 0] 

Putting everything together and doing some simplification, the short distance potential to the two-loop order 
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It is interesting to note that the leading log terms to second order would sum to 

I 
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(A7) 



1 - (33 - 2rc/)(]n(>r) + 7)a s /67r ' 

and that the sub-leading log term can be written 

(31 - 10n//3)as/127r 
1 - (33 - 2n/)(ln(^r) +7)Ss/3tt ' 

As expected, as fir gets small the log term becomes positive and increasingly large thereby softening the potential. 
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Appendix B: Calculational Details 



In the determining the trial wave functions, we use a variational approach as described in |l|. The trial function 
is given by 

n 

= E c * ( r / R ) k+£ e ~ r/R yrM ■ ( B1 ) 

fe=0 

In the determining of the matrix elements using this trial function and the potential given in Eq. (|A6[) . it is worth 
noting that a special type of polygamma function had to be used. The second order logarithmic term of the 
potential when integrated results in an integral of the form 

poo 

/ t z - 1 e- t ln 2 (t)dt = r{z)(¥°\z) 2 + ^ 1 \z)), (B2) 
Jo 

where ^( n \z) represents a polygamma function of order n. In our case we encountered both the digamma and 
trigamma functions. 
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